Transient currents in a molecular photo-diode 
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Light- induced charge transmission through a molecular junction (molecular diode) is studied in the 
framework of a HOMO-LUMO model and in using a kinetic description. Expressions are presented 
for the sequential (hopping) and direct (tunneling) transient current components together with 
kinetic equations governing the time-dependent populations of the neutral and charged molecular 
states which participate in the current formation. Resonant and off-resonant charge transmission 
processes are analyzed in detail. It is demonstrated that the transient currents are associated with a 
molecular charging process which is initiated by photo excitation of the molecule. If the coupling of 
the molecule to the electrodes is strongly asymmetric the transient currents can significantly exceed 
the steady state current. 
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I. INTRODUCTION 

The use of molecular nanostructures as 
diodes, transistors, switches, etc. is considered 
as one possible way towards a further minia- 
turization of integrated circuits. Although pio- 
neering ideas in this direction have been formu- 
lated more than 30 years ago [l|, the detec- 
tion of current- voltage characteristics of a sin- 
gle molecule became only possible within the 
last 15 years (for an overview see [3-10]). Up 
to date research still mainly focuses on an un- 
derstanding of charge transmission in the junc- 
tion "electrodel-molecule-electrode2" (l-M-2 
system) where a single molecule exhibits itself 
as an electron/hole transmitter. It has been 
shown that at definite conditions the molecule 
is able to operate as a molecular diode. For in- 
stance, during a coherent (elastic) electron tun- 
neling in the biased l-M-2 system, diode prop- 
erties of the molecule appear only in the pres- 
ence of a voltage drop across the molecule. This 
conclusion is valid even at different contacts of 
the molecule with the electrodes 111. But, if 



an electron transmission is associated with in- 
coherent electron transfer processes (inelastic 
tunneling or/and hopping), the molecular diode 
can originate from an unequal coupling of the 
molecule to the electrodes. Just such a situ- 
ation is considered in the present paper. We 
show that a rectification effect can be observed 
even in the unbiased l-M-2 system where the 
driving force of the electron transfer process is 
caused by a photo excitation of the molecule. 

Recent research addressed the use of or- 
ganic molecules in molecular photo devices like 



photo-diodes, photo-resistors, optical switches, 
and photo-amplifiers [HI, [l3j]. For example, a 
light-controlled conductance switch based on a 
photochromic molecule has been demonstrated 
[14] . Moreover, single molecule luminescence 
caused by the current through a molecular junc- 
tion could be detected [l5T - tl8j . 

Theoretical estimates on the light-induced 
current and current-induced light emission one 
can find in [l9l - [25| . It has been shown that a 
rfc-current can be induced by an external ac- 
field either due to a considerable difference be- 
tween the electronic charge distributions within 
the molecular orbitals (MOs), or if the ampli- 
tude of the electric field along one direction is 
larger than in the opposite direction. The lat- 
ter effect can be originated by a mixin g of two 
laser pulses with frequencies u) and 2w [2ll . [26l — 
l29j . The generation of a dc-current can be 
also achieved by an asymmetric distribution of 
molecular energy levels caused by environmen- 
tal fluctuations (such an asymmetry may in- 
duce a ratchet current 30]). Besides the for- 
mation of a steady-state current due to an op- 
tical excitation of the junction in the absence of 
an applied voltage, a light-induced suppression 
of a current in the presence of a n applie d volt- 
age has been suggested as well [33, |32| ■ Note 
also the work on a light-induced removal of the 
Franck-Condon blockade in a single-electron in- 
elastic charge transmission [33| . 

While the examples mentioned above focus 
on steady-state properties of the junction, also 
the formation of transient currents generated 
just after an alteration of the applied voltage or 
by changing optical excitation attracted recent 
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interest. The computations demonstrated that 
the transient current in a molecular diode ap- 
pearing just after a sudden voltage switch-on or 
switch-off can significantly exceed the steady- 
state value [H], |38[ . Such a behavior is caused 
by electron transfer processes which are respon- 
sible for charging or discharge of the molecule 
and which are fast compared to the processes 
that establish the steady state current. 

It is the objective of the present work to 
study the time-dependent behavior of the tran- 
sient current across a molecular junction in the 
absence of an applied voltage. In doing so we 
focus on the transient current formed just after 
a fast switch-on of a cw-optical excitation. Our 
analysis allows to clarify the physical mecha- 
nisms which are responsible for the fast and 
the slow kinetic phases of charge transmission 
through the molecular junctions. 

The paper is organized as follows. General 
expressions for the sequential (hopping) and di- 
rect (tunnel) current components in a molecular 
junction are given in the next section along with 
the kinetic equations for the molecular state 
populations and respective transfer rates. In 
sec. Ill, a HOMO-LUMO description of the 
molecule is introduced to derive concrete ex- 
pressions for contact as well as inelastic tunnel 
rates. Expressions for the transient photocur- 
rent are presented in sec. IV. In sec. V, the 
results related to the off-resonant and the reso- 
nant regime of current formation are discussed 
in details. Some concluding remarks are pre- 
sented in Sec. VI. 



II. BASIC EQUATIONS 



a. Electron creation and annihilation operators 
are denoted by a^ ka and a r k<r, respectively. The 
expression 

H m = E M{N) \M(N))(M(N)\ (3) 

defines the Hamiltonian of the molecule, where 
Em{n) denotes the energy of the molecule in 
state \M(N)). The quantum number M labels 
the actual electronic, vibrational, and spin state 
of the molecule; N denotes the number of elec- 
trons in the molecule. The third term in eq. 
([T]) reads 

H' = J2 \ Vm 

'(N+l);rkaM(N) 

rkcr N,MM' 

x \M'(N + l))(M(N)\a rku + h.c] . (4) 

It describes the molecule-electrode interaction 
with the matrix element Vm 1 (N+i);rkaM(N) = 
(M'(N + l)\V tr \r\s.aM(N)) characterizing the 
electron exchange (V* r is the electron transfer 
operator). The interaction of the molecule with 
an external cw-field is written in the standard 
form 

H f (t) = -E(t) d M >M\M'(N))(M(N)\ 

MM'N 

(5) 

where E(i) is the electric component of the pe- 
riodic field and &M 1 (n)m(n) is the transition 
dipole matrix element between different states 
of the molecule. 



A. Hamiltonian 

We introduce a model of the l-M-2 molecular 
junction formed by two nonmagnetic electrodes 
which arc weakly coupled to the (nonmagnetic) 
molecule. The related Hamiltonian of the sys- 
tem can be written as 

H = H e + H m + H' + H f (t) . (1) 

The first term describes the Hamiltonian of the 
ideal electrodes, 

rka" 

where E r ^ denotes the energy of a conduc- 
tion band electron (with wave vector k) of the 
r(= l,2)th electrode. For nonmagnetic elec- 
trodes and in the absence of a magnetic field 
this energy does not depend on the electron spin 



B. Sequential and direct components of 
an electron current 

The current across the electrode r is given by 
I r (t) = e(5 rtl -8 r>2 )N r (t) (6) 

where e = — |e| is the electron charge, N r (t) = 
Skcr -P(»~kcr;£) denotes the electron flow from 
the rth electrode, and P(rka;t) is the popu- 
lation of the single-electron band state. For 
stationary charge transmission the number of 
electrons leaving one of the electrodes is identi- 
cal with the number of electrons arriving at the 
other electrode, i.e. Ni{t) — —N2(t) — const. 
In the nonstationary regime, however, Ni(t) 
and N2(t) may be quite different from each 
other so that h{t) f I 2 {t) (cf. refs. (HHI). 

Nonequilibrium density matrix (NDM) the- 
ory [35l437| is quite suitable to achieve a unified 
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description of elastic (coherent) as well as in- 
elastic (hopping and incoherent) charge trans- 
mission in the molecular junctions. Such de- 
scription allows one to express the transfer rates 
characterizing the noted transmission via the 
set of molecule-electrode couplings and trans- 
mission gaps. In refs. (38Ui^ |. the NDM theory 
has been used to derive kinetic equations for the 
single electron populations P(rker;i) and the 
molecular populations P(M(N);t). Just these 
equations determine the evolution of the cur- 
rent components in the molecular junctions. 

In the presence of an external cu;-field, the 
derivation procedure becomes more compli- 
cated. If, however, the interaction, eq. ©, acts 
as a perturbation only, the calculation of elec- 
tron transfer rates associated with the interac- 
tion, eq. ((H) , can be carried out by ignoring the 
molecule-field interaction. The condition that 
permits one to consider the interaction ([5]) as a 
perturbation, reduces to the inequality 

u 2 > |Ed Afw | 2 A 2 (7) 

where |E| is the amplitude of the cw-field. Ow- 
ing to the condition (|7|), only single photon 
transitions with frequency w = (l/h)\E^j^ N j — 
Em'(n)\ will support charge transfer processes 
in the l-M-2 device. The derivation of kinetic 
equations for the populations P(rka; t) and 
P(M(N); t) remains identical with that already 
presented in [38l442l |. Therefore, we do not re- 
peat the derivation here. We only mention that 
for the considered weak molecule-electrode cou- 
pling a unified description of charge transmis- 
sion is achieved by using the transition operator 
f = H' + H'G(E)H' (note that the electrodes 
stay in equilibrium). The matrix elements 
(a\T\b) determine the transitions between the 
states b and a on the energy shell E = E a = Eb 
[431 ] . where the E a and Ef, are energies refer- 
ring to the Hamiltonian Ho = H e + H m . The 
Green's operator G(E) = (H + H ' + iO+y 1 is 
defined by the Hamiltonian of the whole l-M-2 
system in the absence of molecule-field inter- 
action. The first term of T is responsible for 
a single electron hopping between the molecule 
and the attached electrodes while the second 
term results in a direct one-step electron tran- 
sition between the electrodes. Besides, the op- 
erator H'G(E)H' is responsible for a specific 
electron-pair transition between the molecule 
and the electrodes. Resp ective transfer rates 
are presented in refs. [38l l42j. The mechanism 
of electron-pair transitions has been applied 
earlier to explain the nonlinear electron trans- 
port through a single- level quantum dot (44J. 
For such a transport the repulsion between the 
transferred electrons in the dot is compensated 



by the voltage bias. In the present paper, a 
light-induced charge transmission is considered 
in an unbiased molecular junction. Therefore, 
the energies of twofold charged molecular states 
are arranged high enough to only give a neg- 
ligible contribution to the current. It means 
that the pair-electron transfer processes be- 
come unimportant and, thus, our study is lim- 
ited to single electron transmission processes. 
As a result, the current through the rth elec- 
trode has two components 

Ir(t)=4ll(t)+I dir (t). (8) 

The sequential component, 

N,MM> 



-X ( M\ N ^ MI{N ^)P(M(N);t) (9) 

is defined by single electron jumps through the 
contact region between the electrode surface 
and the molecule. Respective hopping trans- 
fer rates can be referred to the contact for- 
ward (electrode-molecule) and contact back- 
ward (molecule-electrode) rates which read 

(r) _ 27T ^ i T/ i2 

X M (N)^M'(N+l) ~ 1^ \ V M'(N+l);r-kaM(N)\ 

kcr 

x / r (-E r k) S[E r k + Em(n) — ^M'(iV-f-i)] (10) 
and 

X-M(N)^M'(N-l) = \ V M>(N-l)rkcr;M(N)\ 2 
lea 

x [1 - fr(E rk )] 5[E r k + - ^m(jv)] ■ 

(11) 

[In eqs. (HUD and {TTJ, fr(E rk ) = 
{ exp [(Erk — fi r )/ksT] + is the Fermi dis- 
tribution function with [i r being the chemical 
potential for the rth electrode.] Contact for- 
ward and backward rates are responsible, re- 
spectively, for reduction and oxidation of the 
molecule by the rth electrode. In mesoscopic 
physics, a similar type of hopping processes is 
classified as an electron tunneling between the 
lead and the dot [45| . 
The current component 

I dir (t) = \e\ Yl S^,(N)P(M(N);t) 

N.MM' 

(12) 
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is formed by an interelectrode electron transfer 
at which the molecule mediates a charge trans- 
mission without alteration of its charge. Such 
process is defined by the electron flows 

Sffl'M = QlM(N)^2M'(N)—Q2M(N)^>lM'(N) 

where the transfer rates 

QrM(N)-yr'M'(N) 



kcr, k'cr' 



x\(M'(N)r'k'a'\H'G(E)H'\rkaM(N)) 



X S[E rk + E M (N) ~ #r'k' - E W (N)] (14) 

characterize a distant electron transmission 
from the kcr band states of the rth electrode 
to the kV band states of the r'th electrode. 
Such transmission appears as a direct single- 
step elastic (at M'(N) = M(N)) or inelastic 
(at M'(N) ^ M (N)) interelectrode electron 
tunneling. Since the operator H' is respon- 
sible for transitions accompanied by an alter- 
ation of molecular charge, the mediation of the 
tunneling transmission occurs via the forma- 
tion of intermediate molecular states M(N + 1) 
and M(N — 1) which differ from the initial, 
M(N) and final, M'(N) charge states. In the 
contrast to the sequential (hopping) transmis- 
sion where similar states are really populated, 
the noted intermediate states are not populated 
and only acts as virtual states. In mesoscopic 
physics, such type of transmission refers to co- 
tunneling 45]. In the respective terminology 
the direct current component, eq. (1T3|) . results 
as a contribution of partial currents associated 
with different co-tunneling channels. The real- 
ization of a particular channel is controlled by 
the probability P(M(N);t) to find a molecule 
in the M(iV)th stay. [Examples of electron 
transmission along the channel pathways that 
include the empty and occupied MOs can be 
found in [H, S3].] Thus, the direct tunneling 
can be referred to as co-tunneling which is con- 
trolled by kinetic charging and recharge of the 
molecule (via electron jumps through the con- 
tact region). This circumstance has been al- 
ready noted in 48] . In the presence of the cw- 
ficld, an additional control occurs through the 
population of the excited molecular state. 



C. Kinetic equations for the molecular 
populations 

It follows from eqs. ([9]) and (fT2)) that each 
charge transmission route (sequential or direct) 
includes electron transfer channels related to 
the molecular states M(N). The contribution 
of the M(iV)th channel to the route is weighted 
by the molecular population P(M(N); t), which 
satisfies the normalization condition 

Y / P(M(N);t) = l. (15) 

NM 

Follo wing the derivation procedure presented in 
refs. (Hi [49[ and bearing in mind the fact 
that the interactions @ and (JSJ are considered 
as perturbations, we can see that evolution of 
the P(M(N);t) is determined by the balance 
like kinetic equation 

P(M(N);t) = -J2 [^M(N)^M'(N')P(M(N);t) 

M'N' 

JC M , m ^ MW P(M'(N');t)} . (16) 
The transfer rate 

£m(n)-*M'(n') — (5n',n+i + 6n',n-i) 



X Xm(N)->M>(N>) +$N : N> \ K M(N)^M'(N) 

rM(N)->r'M'(N')] (17) 

r 

specifies the transition from the state M(N) to 
the state M'(N') in the molecule. Such transi- 
tion is caused by the molecule-electrode interac- 
tion Q through the contact and distant trans- 
fer rates (eqs. flUB CU}, ([T3|l) as well as by the 
molecule-field interaction Respective rates 
of optical excitation and de-excitation are 

K-W - 2n IFM I 2 

M(N)M' (N) ~ ~ft \^ a M'(N)M{N)\ 

x [Lm'(n)m(n)(^>) + Lm'(n)m(n)(~^)} ■ (18) 
We introduced Lmi(n)m(N){w) = 

(1/27T) (k M (N) + KM'{N)){[fiM - (E M '(N) ~ 

E M (N))] 2 + (km(n) + K Af(Ao) 2 /4}~\ where 
K M(Af)/2 denotes the molecular level broaden- 
ing caused by electron-phonon interaction as 
well as interaction of the molecule with the 
electrodes (for more details see [49[). 



5 




FIG. 1: Possible position of the frontier MOs in 
the molecule with two terminal sites I and II. In- 
tersite coupling transforms the site MOs into ex- 
tended HOMO, HOMO-1 and LUMO, LUMO+1. 
Spots indicate the main location of electron density 
within the extended MOs. 



III. CHARGE TRANSFER PROCESSES 
IN THE HOMO-LUMO MODEL 

Next, the hopping rates, eqs. ([TU]) and (ITTj) 
as well as the distant transfer rate, eq. (ITU) , all 
determining the net electron flow through the 
junction, are further specified along the rate, 
eq. (fT5|) characterizing the efficiency of exci- 
tation and de-excitation of the molecule. We 
use a model of the l-M-2 system where only 
the highest occupied and the lowest unoccupied 
molecular orbitals (HOMO (H) and LUMO 
(L), respectively) are considered. The HOMO- 
LUMO model is suitable to study charge trans- 
mission in the molecular junctions. As an ex- 
ample, note the pioneer work of Aviram and 
Ratner [50| where the mechanism of current for- 
mation includes a participation of HOMO and 
LUMO levels belonging to the donor and ac- 
ceptor sites of the molecule. Recently, a similar 
model (with chromophoric donor and acceptor 
sites) has been used for the description of tran- 
sient dynamics in a molecular junctions [511 ] . 
In this model, a transient electronic current is 
formed due to an optical excitation associated 
with the HOMO-LUMO transition in the donor 
site. 

In the present paper, we use a model where 
extended HOMOs (LUMOs) are formed from 
HOMOs (LUMOs) belonging to the terminal 
molecular sites I and II coupled to one another 
by interior bridging groups. Let the HOMO(ra) 
and the LUMO(n) refer to the terminal site 
n(=I,II). Following from the coupling between 
the sites, the extended HOMO and HOMO -1 
(LUMO and LUMO +1) represent a mixture of 



the HOMO (I) and HOMO (II) (LUMO (I) and 
LUMO (II)). If the intersite coupling does not 
strongly modify the electron distribution across 
the molecule, the maxima of electron density 
in the HOMO, HOMO -1, LUMO, and LUMO 
+1 correspond to electron densities located in 
the vicinity of the respective sites, cf. Fig. 
[0 Therefore, the coupling of the HOMO to 
electrode 1 is assumed to be much stronger 
than the similar coupling to electrode 2. The 
opposite case is valid for the coupling of the 
LUMO to the same electrodes. This configu- 
ration as represented in Fig. [T] can be realized 
if, for instance, the HOMO(I)/LUMO(I) and 
HOMO (II) /LUMO (II) refer to the 7r-electrons 
of aromatic groups coupled to each other by the 
bridging cr-bonds (to avoid a noticeable mixture 
between the 7r-electrons belonging the sites I 
and II). If the energy hco of the external cw- 
field coincides with the energy of the optical 
HOMO-LUMO transition, then the formation 
of the photocurrent can be mainly associated 
with two frontier MOs (HOMO and LUMO). 
In this case, the rates of optical excitation and 
de-excitation are determined by eq. ([3]). 

For the subsequent analysis we assume that 
the Coulomb interaction between excess elec- 
trons (or holes) occupying the molecule in the 
course of charge transfer, is so large that the 
molecule can only stay in its neutral ground (or 
excited) state, in its oxidized state and in its re- 
duced state. These states are denoted as M$ = 
M(N G ), M* = M'(N G ), M+ = M(Nq - 1) 
and M_ = M(N G + 1). Here, N G is the num- 
ber of electrons if the molecule is in its neutral 
state. If the maxima of electron location at the 
HOMO and LUMO are in the vicinity of the 
spaced sites I and II, cf. Fig. [TJ one can sup- 
pose that the exchange interaction between the 
unpaired electrons occupying the HOMO and 
the LUMO becomes small. This allows one to 
ignore the exchange splitting between the sin- 
glet, M*(S) and triplet, M*(Tm),(m = 0,±1) 
states of the excited molecule. Accordingly, 
the electron spin projections can be taken as 
good quantum numbers. Therefore, the four- 
fold degenerated excited state M„ can be char- 
acterized either by molecular spin states (M* = 
M* (S), (Tm)) or by spin projections an and 
(7£ of unpaired electrons occupying the frontier 
MOs (M* = M*{a H ,a L )). At a negligible ex- 
change interaction, both sets of spin quantum 
numbers lead to identical results. Moreover, 
the states M+ = M+(a H ) and M_ = M_(cr L ) 
are twofold degenerated. 

To specify the energies E M ^) entering the 
molecular Hamiltonian and the matrix ele- 
ments in the interaction expression §Q we intro- 
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duce the following notation of the Hamiltonian 

j " 

+ \ E J2 u ^'4-' c r-') c t c ^ ( 19 ) 

and of the electron transfer coupling (cf. flU 

j rker 

(20) 

In eq. (H?)) the ej are the energies of an electron 
occupying the j(— H, L)th MO. The strength of 
the Coulomb interaction between two electrons 
is defined by Uj if both electrons occupy the jth 
MO. If the electrons belong to different MOs 
Coulomb interaction is measured by Ujji . The 
operators Cj a and Cj a create or annihilate an 
electron in the molecule, and /3j>k characterizes 
the coupling of the jth MO to the rkth band 
state of the electrode. 

According to the Hamiltonian, eq. (|T9|). 
the molecular energies Emin) — E<x: i a = 
0, *, +, — ) follow as: 

E = 2e H + U H , E„ = e H + e L + U LH , 

E-=26 H + e L + U H + 2U LH , E+=e H . (21) 

Here, tu and €l are the energies of an 
electron occupying the frontier MOs while 
Uh and Uhl are the Coulomb parameters. 
Concerning the matrix elements entering eq. 
(|4|). all of them are expressed by the cou- 
plings pHrk or 0Lrk- We have, for exam- 
ple, (M \V t r\M + (a H )rka) = l3 H rk$-a.a H and 
(M rka\V tr \M_(a L )) = /?£ rk <W 

A. Contact rate constants 

Noting the structure of the transition ma- 
trix elements, the so-called wide band approx- 
imation [55| enables one to express the hop- 
ping transfer rates (fT0|) and (jTTJ) by contact 

(r) 

rate constants K ',. For instance, we get 

XM, { a H ,a' L )^KU(a L ) = ^ L ,>n. K i-- According 
to the used HOMO-LUMO model the forward 
contact rate constants takes the form 

K^c(l/h)r^N(AE-o), 
Ki r ±~(l/H)rPN(AE_*), 



K^ ^(l/h)rPN(AE 0+ ), 

K { ;1 ~ {l/h)T { [ ] N{AE*+) . (22) 
The quantities 

k 

characterize electron hopping between the jth 
MO and the rth electrode (cf. Fig. [2|), whereas 
the distribution function 

N(AE a , a ) = [exp (AE^Jk B T) + lp 1 . (24) 

determines the influence of temperature on the 
hopping processes via the transmission gaps 

A£'_|_o(*) = (E+ + Ep) - i?o(*) (25) 

and 

A£L oW =E_- (E 0H + E F ) . (26) 

Backward contact rate constants which char- 
acterize the transition of an electron from the 
molecule to the rth electrode are connected 
with the forward ones, eq. (|2"2"|) . by the rela- 
tion 

K { :l = K$ a exp (~AE a , a /k B T) . (27) 

The physical meaning of the transmission gaps 
can be easily deduced from their definition. 
Since E- — E and E + — Eq are the elec- 
tron charging and electron discharging ener- 
gies (with respect to the molecule being in its 
ground neutral state), respectively, the inequal- 
ities E- — Eq > Ep and E — E + < Ep have to 
be fulfilled in the unbiased l-M-2 system with 
identical electrodes (cf. the upper panel in Fig. 
[3] where Hi — ^2 — Ep). Therefore, the gaps 
A_E_o and AE +n are both positive. When the 
molecule is in the excited state, then respec- 
tive charging and discharging energies, £L —E* 
and E + — E*, can be higher or lower than the 
Fermi level and, thus, transmission gaps A_E_* 
and AE+* can become positive or negative. 
[One possible case with E_ — E* > Ep and 
E* — E + < Ep, is presented at the upper panel 
in Fig. O] 

An additional interpretation of the transmis- 
sion gaps follows from a comparison of elec- 
tron energies belonging the whole l-M-2 system 
[H,[3^|. Let £ e be the energy of electrons in the 
electrodes. In the case of a neutral molecule the 
energy of the whole system is E{\Mqu\2) = 
Eqm + £ e . During a charge transmission pro- 
cess the number of electrons in the system is 
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FIG. 2: HOMO-LUMO scheme related to the elec- 
tron transfer through the l-M-2 molecular junc- 

(r) 

tion. The width parameters Tj characterize the 
efficiency of contact electron jumps as well as of the 
direct (tunneling) electron transfer. The sequential 



current components lill(t) and Istq{t) can differ 
from each other if the junction is transient regime. 



(2), 



E E 4 
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conserved. Therefore, the energies of the sys- 
tem with the oxidized and reduced molecule 
are, respectively, E(l~M + 2) = £(lM + 2") = 
E + +£ e + E F and£ , (l+M_2) = £(1M_2+) = 
E_ + £ e - E F . Therefore, the gaps (J25) 
and ([2"B)) correspond to the difference between 
the above noted energies, i.e. Ai? +0 (*) — 
E(l~M+2) - E(1M 0( *)2) = E(lM+2-) - 
E{lM m 2) and AE_ 0M = £7(1+M_2) - 
£?(1M 0W 2) = £(1M_2+) - £7(1M 0W 2) (cf. 
Fig. [3] lower panel). Such interpretation of the 
transmission gaps is quite suitable for the anal- 
ysis of the transmission processes in the molec- 
ular junctions. 

The sign of the transmission gap defines the 
electron transfer along a given transmission 
channel. For instance, if A_E + * is positive, then 
the transition M + —> M* caused by an electron 
injected into the molecule, requires a thermal 
activation i.e. it proceeds in an off-resonant 
regime. If AE + * < 0, however, the M + — > A/* 
transition does not require any thermal acti- 
vation and, thus, becomes practically indepen- 
dent on the absolute value of A_E + ». Conse- 
quently, the electron hopping takes place in a 
resonant regime. 



B. Inelastic tunnel rate constant 

Eq. (|14p for the direct (tunnel) transfer rate 
indicates that in the absence of an applied volt- 
age any elastic electron tunneling between iden- 
tical electrodes disappears. We study in the 
following an inelastic tunneling event which is 
accompanied by the intramolecular transition 
M* — > Mq. To derive a respective rate expres- 
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FIG. 3: Charging energies E -E+, E* -E+, E-E 
and E- — E* (upper panel) and transmission gaps 
related to the charged molecular states M+ and M- 
(lower panel). 



sion we first consider the effect of the molecule- 
electrode coupling. Since it is not too strong, 
its presence may be accounted for by a shift 
AE M i jv) of the molecular energies as well as by 
a level broadening 

rA/(AT)/2 = {\yM(N)rka:,M'(N+l)\ 2 
r M' ko- 

+ \^M(N);M'{N-l)rkcr\ 2 S[EM(N)— ^M'fJV-l) — -^rk]} 

(28) 
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Accordingly, the Green's operator, determin- 
ing the general transition amplitude is defined 
by these shifted and broadened molecular ener- 
gies. Due to the weak molecule-electrode cou- 
pling one can omit the energy shift A.Em(n) ■ 
Taking the Hamiltonian, eq. ()19[) . then, in the 
framework of the HOMO-LUMO model, the 
unperturbed molecular energies, eq. (j2"Tj) arc 
expressed via the single-electron energies ch 
and £i as well as via the Coulomb parameters 
Uh and Ulh- The broadenings, eq. (|28| are 
defined by the single-electron level broadenings 
Tj/2. The quantities 



(r) 
j 



U = H,L). 



(29) 



are obtained as the sum of the width param- 
eters r^- r ', eq. P5|) . Thus, the described for- 
mulation of the Green's operator G{E) and the 
introduction of the width parameters results in 
the following rate expression 



!rAf,((Tt,(7H)->r'Jlfo 



r (r') r (r) 
L 1 H 



1 



r (r') r (r) 



(30) 



used 

,(r) 



r 4 



(r) 



2l£>) 



and 



Here 

r_ = J2r ( t l ] + 2r ff)- According to eq. fl30 
the regime of inelastic tunneling transfer is gov- 
erned by the quantities 



arc tn (2AE a > a /T a i 



(31) 



For the weak molecule-electrode coupling un- 
der consideration the width parameters do 
not exceed 10~ 3 eV. Accordingly, we have 
\AE a > a /T a >\ 3> 1, and one can use the 
asymptotic form <p a <^ a ~ (7r/2)(signA-E a ' a ) — 
(T a r /2AE a / a ). It follows the particular rela- 
tion (p+^o — ~ (7r/2)(l — signA_E + ») + 
[(r+/2AP+*) - (T + /2AE +0 )} » (tt/2)[(1 - 
signAP+*) + (r+/7rAE + *)]. Taking AE+* > 
then ip+^o — « (r + /7rAi? + ») C 1. In 
the contrary case A_E + » < we find ^+^o 
.,: ^ ?a 7r what is much larger then the respec- 
tive expression deduced for AE +tf > 0. Note 
also that at the resonant regime of tunnel elec- 
tron transmission the ip^ >o — f-i >* is inde- 
pendent of the actual value of the transmission 
gap. 



C. Rate equations for integral molecular 
populations 

The kinetics in the junction are dominated 
by sequential processes which are character- 

(r) 

ized by contact rate constants K^a" But, * n e 
direct inelastic electron tunneling responsible 
for the transitions between the excited and the 
ground states of the neutral molecule, is also 
important. For instance, the distant rate con- 
stants describe the nonradiative decay of 
the excited molecule. The kinetic scheme de- 
picted in Fig. 2] illustrates the possible tran- 
sition processes in the junction including the 
cw-optical excitation of the molecule. All rates 
indicated in Fig. @] characterize the transi- 
tions with the participation of the degenerated 
molecular states M*, M + and M_. Therefore, 
it is convenient to introduce integral molecular 
populations 



P(*;t)= P(M*(*H,<T L );t), 



P(-,t) = J2P(M-(vL);t),. (32) 



[Note also that P(*;t) = P*(S;t) + 
J2m=o ±1 P*(T m )'i t)-} The quantities intro- 
duced in eq. (|32|) along with the population 
P(0; t) = P(Mq; t) obey the normalization con- 
dition 



a— 0,*. + , — 



(33) 



which corresponds to eq. (fT5|) . Based on the in- 
troduction of integral populations, the general 
kinetic equations (fTo| reduce to the following 
set of rate equations 

P(0; t) = ~k f P(0; t) + K +0 P(+; t) 



+K- P(—,t) + k d P(*;t) , 



P(+;t) = -{K +0 +2K +m )P(+;t)+K^+P(*;t) 



P(-; t) = -(#_„ + 2K_*)P{-; 1) + K^P*(t) , 



P(*; t) = -(#*+ + + k d )P(*; t) 
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FIG. 4: (color online) Kinetic scheme of the trans- 
fer processes occurring in the molecular junction in 
the absence of an applied voltage (for further de- 
tails see text). 

+2K +m P(+;t) + 2K-J>(-;t) + kfP (t) . (34) 
Here, we introduced the recharge transfer rates 

Kaa^KM+KW (35) 

which are expressed by the sum of contact 
rate constants. The rate constant kf = 2£q{ , 
cq. (fl~8f characterizes the optical transition be- 
tween molecular singlet states M(N) — Mo and 
M'(N) = M»(S). Accordingly, the overall de- 
cay rate from the four-fold degenerated excited 
state follows as 

k d = kf/A + Q^ . (36) 

We defined 

Q*o = Qio 2) + Qil 1] (37) 

as the component caused by the coupling of the 
molecule to the electrodes. The concrete ex- 
pression for Q* can be deduced from eq. (|30|). 
It reads 

fto^^f + rMV (38) 

with 

R = r; 1 [(1 - signA£+*) + (T + /nAE + *)} 
+ n 1 [(1 - signAEU) + (T_/7rA£L„)] . (39) 

IV. THE PHOTOCURRENT 

Already in the absence of an applied voltage 
a photocurrent, eq. © has to be expected. Ac- 
cording to the used HOMO-LUMO model and 
by noting the eqs. ([TO ]) . ([TT ]> . ([22 ]) and ([27 ]) for 



the sequential current component, eq. ([9]), we 
find 

iilKt) = (<Jr,i -5 r , 2 )I nh{[Ki r 2P(*;t) 
+(K^+2K^)P{+;t)]-[K^P(*;t) 

+ (K^+2K^)P(-,t)}}, (40) 

where Iq = \e\/ivh x leV ~ 77.6 /xA is the cur- 
rent unit [56|. It follows from eq. (l40l) that the 
charge transmission along the sequential route 
is determined by hopping (contact) rate con- 
stants ([22)1 and ([27)1 . The time-dependent be- 
havior of this current component is determined 
by the molecular state populations P(+;t), 
P(—;t) and P(*;t). The direct component of 
the current is formed by the inelastic tunnel 
electron transmission along the channel related 
to the excited molecule. The expression for the 
direct component follows from eqs. ([12")) , (|T5|) . 
and ([30|) and takes the form 

(t) = I TThS* P(*;t). (41) 

Here, we introduced 

o _ n (^) ,q(21) 

-^(r^rf-rM 1 ^ (42) 

what represents the net tunnel electron flow (R 
has been introduced in eq. (|5§)l ). Eqs. (l4Tt 
and (|42|) show that, the sign of the direct cur- 
rent component (in the HOMO-LUMO model) 
is determined by the sign of T^T^ - T^r^. 
Moreover, the time-dependent behavior of the 
direct current component is determined by the 
population P(*; t). 

The scheme of the electron transfer routes, 
as displayed in Fig. [SJ offers the opportunity 
to analyze further details of the current for- 
mation. Charge transitions are represented by 
12 — > l + 2~. According to Fig. [5] a possi- 
ble current formation results from the decay of 
the excited molecular state. This is possible 
via the sequential route (including the forma- 
tion of the charged molecular states M + and 
M_) as well as via the direct route M* — > Mq. 
Both routes include two (left and right) trans- 
mission channels. The left sequential chan- 
nel, 1M*2 ^ 1M+2" -> 1+M 2~, proceeds 
across the charged molecular state M + . This 
state is formed by an electron hopping from 
the LUMO to the electrode 2 (with contact 
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FIG. 5: Transfer routes of the light-induced in- 
terelectrode 12 — >• l + 2~ electron transmission in- 
cluding the participation of the charged molecular 
states M+ and M_ . Within the transmission along 
the sequential route, the charged states are popu- 
lated while the same states participate in a virtual 
form only if the transfer proceeds along the tunnel 
route. 
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FIG. 6: Feasible energy levels of the molecular 
junctions in the absence of an applied voltage. 
Eo, E*, E+ and E- are the molecular energies and 
Ef denotes the electrode Fermi energy. 



change of the charged molecular states popu- 
lations (these states only participate as virtual 
states). Thus, a direct electron transfer consti- 
tutes an inelastic tunneling event of an electron 
between the electrodes. The related distant 
transfer rate Ql is defined in eq. ([50)1 . Since 
a reverse route 1Mq2 4— l~A/*2 + is formed in 
a similar way and is characterized by the rate 

(21) 

Ql Q , the direct current component is propor- 
tional to the net electron flow S*q, eq. P^j). 



RESULTS AND DISCUSSION 



rate constant K#l) after which another elec- 
tron hops from the electrode 1 to the HOMO 
(with contact rate constant K+q). In sum- 
mary, we have the transition 12 — > l + 2~. 
This transition is also achieved by an electron 
transmission along the right sequential chan- 
nel 1M*2 <=► 1+M_2 -> l + Mo2" which in- 
cludes the charged molecular state M_. The 
M_ state is formed by an electron hopping from 
electrode 1 to the HOMO. The transfer of an 
electron from the LUMO to electrode 2 returns 
the molecule to its neutral ground state Mq. 
The respective sequential charge transfer steps 
are characterized by the contact hopping rates 
K±z and K_\. It is important to underline 
that during the electron transmission along the 
sequential route, the current formation is ac- 
companied by a molecular recharging, i.e by a 
population of the intermediate charged molec- 
ular states M + and M_. 

The second type of transfer route shown 
Fig. [5j refers to the direct (distant) interelec- 
trode electron transfer 12 — > l + 2~ which is 
accompanied by the M* — > Mq transition in 
the molecule. In contrast to the sequential 
route, electron transmission along the direct 
route 1M*2 — > 1 + Mq2~ does not result in a 



The eqs. ©, (gDJl, and (gj) allow one to 
describe the time-dependent evolution of the 
photocurrent in the molecular junction start- 
ing with the switch-on of a cw;-optical excita- 
tion (at t = 0) and extending up to the for- 
mation of a steady-state current I st = I\(t 
T st ) = 22 (i S> T s t) where r st is the character- 
istic time of the steady state formation. The 
kinetic schemes drawn in the Figs. 2] and [5] 
show that a control of the light-induced electron 
transfer is achieved via the transitions between 
electronic states Mo,M*,M + and M_ of the 
molecule. Here, the charged molecular states 
M+ and M- are of particular importance since 
those participate in the transmission channels 
formation related to the sequential and the di- 
rect electron transfer routes. As far as the 
charged states population is determined by the 
relation between forward and backward contact 
rate constants, eqs. (PZ21 and ([2"7]). the position 
of the energy levels in the l-M-2 system prede- 
termines the specific form of the interelectrode 
12 ->• l+2~ and 12+ <- 12 electron trans- 
mission. 

Possible arrangements of the molecular junc- 
tions energy levels are depicted in Fig. [6] For 
the sake of simplicity, the electron energy £ e 
of the electrodes is omitted (cf. Fig. [3J lower 
panel where this energy is presented). If the 
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junction energy with the neutral molecule in 
its excited state is below the energy valid if the 
molecule is in a charged state, then only an off- 
resonant regime of light-induced electron trans- 
mission becomes possible (cases (a) and (b)). 
The cases (c) and (d) correspond to charge 
transmission with a single resonant channel (ei- 
ther the oxidized or the reduced molecule is in- 
volved) . Two further transmission channels are 
realized if the junction energy with the neu- 
tral molecule in its excited state exceeds the 
energies realized for the oxidized or the reduced 
molecule (cases (e) and (f)). 

The further analysis will be based on the eqs. 
©, gO]), (HU) for the current as well as the 
eqs. (|22[) and (|27|) defining the contact rate 
constants. Additionally, the expressions (157)) 
- for the tunnel decay rate Q* and eq. 
(|4"2")) for net electron flow S*o are taken into 
account. The time-dependent evolution of the 
integral molecular populations P(a; t) are de- 
termined by the rate equations (|M|) where the 
recharge rate constants are given by eq. (|35[) . 
Initial conditions for the populations are found 
from a solution of the rate equations (|3"4l if one 
sets P(a;t) = and kf = 0. Since the rela- 
tions Kq_ ~ and Kq + ~ are valid for the 
case of a charge transmission in the absence 
of an applied voltage, it follows P(*; 0) ~ 0, 
P(+;0) ~ 0, P(-;0) ~ and P(0;0) ~ 1. As 
already indicated the hopping and the tunnel 
transition processes can proceed in off-resonant 
or a resonant regime depending on the sign of 
the actual transmission gap, cf. eqs. (|2"5)) and 

Although the calculations have been per- 
formed on the basis of the general expressions 
(131, and (HOD - gU), including eqs. (|30j), (jM]) 
and (f3"6")) - (|39[) . most of the findings will be 
discussed in terms of analytical expressions. 
Those are derived for cases where an electron 
transfer occurs preliminary along a separate 
transmission channel. For the sake of definite- 
ness, we consider a charge transmission pro- 
cess where the energy P_ — Ep is larger than 
P* and E + + Ep (see the cases (a) and (c) in 
Fig. [5|). We also suppose that AP_* is large 
enough to neglect the population of the state 
M _ . This means that the off- resonant and reso- 
nant regimes of light- induced current formation 
involve an electron transfer process predomi- 
nantly across the three molecular states, Mo, 
M* and M + (the left route in Fig. [5J. Thus, 
charge transmission only occurs along the chan- 
nel related to the charged molecular state M + . 
To achieve an analytic description of this trans- 
mission process one has to set ~ 0. This 



leads to the following solution of the eqs. (|34|) : 

P(0;t)~P Q + £ 

kiK2(k 1 - k 2 ) 

x[k 2 (k 1 -\ + -K* + )e- klt -k 1 (k 2 -\ + -K*+)e- k2 

P(*; t) ~ P.+ % [-fe(fci-A+) e~ fclt 

k 1 k 2 (k 1 - k 2 ) 1 

+k 1 (k 2 - A+) e~ fc2t , 

k\k 2 {ki - k 2 ) 
P(-;i)~0. (43) 

The quantities 

P = (K*+K+o + A: d A + )/fcifc 2 

P+ = k f K*+/k 1 k 2 , P_=0, 

P* = fc/A+/fcifc 2 (44) 

are steady state populations and the overall 
transfer rates take the form 

fci,a = (1/2) [a ± \Ja 2 -4b 2 ] . (45) 

Note the abbreviations a — kf + A+ + A*, 
b 2 = kf(\+ + K* + ) + \*K +0 + 2k d K + *, and 
A + = K +0 + 2K + * , A* = K* + + kd Based on 
the eqs. (14*01) - (|4*3"j) one can derive analytic ex- 
pressions for the sequential and direct current 
components. The time evolution of these com- 
ponents is determined by the overall transfer 
rates k\ and k 2 . Obviously, the formation of a 
finite photocurrent in the absence of an applied 
voltage only becomes possible at an asymmet- 
ric coupling of the molecule to the electrodes. 
As it was already noted such asymmetry can 
result from nonidentical electron density at the 
HOMO and the LUMO (cf. Fig. [T}. For the 
following we assume T% ] > ,T { ^ > . 

Since the factor T^T^ -T^T^ becomes pos- 
itive the steady state electron current is also 
a positive (electrons move from electrode 1 to 
electrode 2). 

Next let us consider the current formation 
at a weak molecule-electrode coupling (when 
the width parameters are of the order (10~ 7 — 
10~ 4 ) eV) and at a moderate optical excita- 
tion (kf = 10 10 s _1 ). The excitation energy is 
fku = E*—Eq — 1.6 eV. The various parameters 
are collected in Table 1. 
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Table I 

Parameters of the HOMO-LUMO model 
(AE±* and are given in eV) 
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A. Off-resonant regime of charge 
transmission 

Figs. [7] - [9] demonstrate the transient behav- 
ior of the direct and sequential current com- 
ponents for a positive transmission gap AE+* 
(it determines the efficiency of the M* — > M + 
transition due to thermal activation). 



1. Deep off-resonant regime 

In the case represented in Fig. [3 the gap 
AE + * is assumed to be a rather large so that 
the population of the charged molecular state 
M + is small. As a result, charge transmis- 
sion mainly follows the tunneling route whereas 
thermal activation of the molecular state M + 
(as well as state M_) is suppressed (cases (a) 
and (b) of Fig. |6]). For such a transmission the 
distant current component strongly exceeds the 
sequential one (see the insert). Therefore, the 
total current, eq. (jSJ is associated with the di- 
rect component: 



7 r (i) ~ 7 dir (t) ~ 7^> (1 



-t/Ts 



(46) 



Here, T s t = k~ t = (kf + fed) is the character- 
istic time of the evolution of the current to its 
steady state value 



T dfJ = IoTthS*o (kf/k s t) 



(47) 



The quantities S*o and Q*o (the latter enters 
k d , cf. eq. (HJoTO have been defined in the eq. 

and eq. (pflZj) , respectively, additionally us- 
ing R ~ (1/tt)[(AE + ,)- 1 + (A7U)" 1 ]. 

In the off-resonant regime, the inequality 
kf 3> Q*o is valid with a good accuracy. As 
a result, a dependence of the photocurrent on 
the light intensity (i.e. on the rate kf) is only 
present within the transient behavior whereas 

the steady state value, I st = Jj**' , becomes in- 
dependent on kj (see Fig. 7). We also note 



FIG. 7: Transient current components for an off- 
resonant charge transmission process and at a dif- 
ferent intensities of optical excitation. Solid lines: 
k f = 10 10 s~\ dashed lines: kf = 3 • 10 10 s _1 . 
Insert: sequential current components /ieg(t) (up- 
per curves) and IsVq it) (lower curves) approach the 
common steady state value. 
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FIG. 8: Off-resonant regime of current formation 
at a small transmission gap AE+*. 



that the single-exponential kinetics correctly 
describes the transfer process until the popu- 
lation of the charged molecular state M + (and 
M_) becomes so small that the direct current 
component strongly exceeds the sequential one. 



2. Single- channel off-resonant regime 



The Figs. I51IT01 depict how the current com- 
ponents approach their steady-state values for 
the case of a large and positive gap A£L* and 
a positive but not so large gap AE + *. Here, an 
activation of the hopping process M* — >• M + 
is possible and the contribution of the sequen- 
tial component to the total one becomes sig- 
nificant. As a result, the transient behavior 
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FIG. 9: Off-resonant current formation at low 
temperature with the participation of the charged 
molecular state M+. 



guarantees k\ 3> fc 2 where k\ ~ kf + kd and 
fc 2 ~ (1/fi) [r L (2 -N{AE+.)) + r H ], Due to 
the inequality fcj 3> &2 the fast and the slow 
kinetic periods of the time evolution are well 
determined. This allows one to distinguish be- 
tween the currents lil q (t) and (t) . The fast 
kinetic phase covers a time region of the order of 
fcj -1 and starts just after the switching on of the 
optical excitation. The phase ends at t > Sk^ 1 
from which the time-dependent behavior of the 
current is determined by the slow phase 



j(*)(l- e -*/r-) +J (rjO») e -*/r rt 



(50) 

Here, T s t = fc 2 is the characteristic time the 
current needs to achieve its steady state value. 
The expression 



of the current represents two-exponential kinet- 
ics. For instance, a comparison of Fig. [7] and 
Fig. [8] shows that at a small gap AE +:t (but at 
the same width parameters and temperature), 
an electron transmission along the sequential 
route becomes more effective than the trans- 
mission along the tunnel route. Therefore, the 
total current is not caused by the direct com- 
ponent (as in Fig. [7|) but by the sequential 
one (cf. Fig. [SJ. As a result, the total cur- 
rents lift) ~ Iseq(t) and ft) ~ /ieg(i) do not 
coincide in the transient region. Such a behav- 
ior is originated by an asymmetric hopping of 
electrons between the molecule and the elec- 
trodes. A temperature decrease does not pre- 
dominantly impact the direct component but 
strongly reduces the sequential one, as deduced 
from a comparison of Fig. ® and Fig. [5] 

For a further inspection of the transient cur- 
rent either in the off-resonant or the resonant 
regime we use analytic expressions which fol- 
low from the eqs. (|30]). (l4Tj) and Let us 
start with an analysis of the sequential current 
components 



Is eg ft) Iseq \}- 



1 



k\ - h 



-k 2 t 



(r) 

+ (-i) r Io7rp^ r N(AE +r )(e- 

fcl — K2 

where the quantity 



-k'2t — k±t\ 



n e q = 7 ° 7r ^^ ( - i H L H L ) N 



denotes the steady state sequential current en 
tering lil\(t) and lil\ft). 

The results depicted in the Figs. H] and O 
refer to an electrode-molecule coupling which 



= (-l) r /o<'iV(AB +t ) 



(51) 



gives the current component valid at t <C r st . 
The sign of the iilq 71 ^ is determined by the di- 
rection of electron motion from the photoex- 
cited molecule to the rth electrode. Since in 
the case under consideration, the transmission 
channel is associated with the charged molecu- 
lar state M + , the — »■ M + transition involves 
an electron which leaves the LUMO and is cap- 
tured by either the 1st or the 2nd electrode. In 
the scheme depicted in Fig. \5\ this transition is 
accompanied by an electron hopping from the 
LUMO to electrode 1 and is characterized by 



(i) 
*+ ■ 



the contact rate constant K 

As a quantifier of the transient kinetics we 
introduce the ratio 



n 



(r) — I j( 
sen Is 



' on) 1 ij(st) 



sc.q 



seq 



(52) 



It indicates how strongly the sequential compo- 
nents differ from their steady state value if the 
optical excitation is switched on. Our studies 
show that the difference between these quan- 
tities can become large if the difference be- 
tween the width parameters is large. We illus- 
trate this observation for the case Tjj 3> 1^ 
and T { 1 ] > r^ 1} so that T^T^ 



r g> r a> 



r (i) r (2)^ Tnerefore; if 5 for instance, = 

o.irf then ijW ~ (r^ r) /r^ ) )(2 - N(AE + *)) 



and thus \I. 

Jst) 



(l,on) I 



< lieq whereas Iltq""' > 



(2, on) 



(48) JW , cf. Fig. E 



During the slow kinetic phase the behavior 
of the direct current component is described by 
the expression 



(49) / dir (^/W(l- e -^) + 4»)e-^ 



The two current components, 



lW=I nhS* 



1 - w 2 N{AE+ *i 



(53) 



(54) 
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and 

Id™ )=I o^S*o, (55) 

represent the steady state values, respectively. 
The basic difference in the behavior of the di- 
rect and sequential current components is as 
follows: In the off-resonant regime of elec- 
tron transfer (where N(AE+*) <C 1) the maxi- 
mal value of the direct switch-on current coin- 
cides with its steady state value, 7"^°™^ ~ I d f r . 
Therefore, a slow kinetic phase is not observed 
for the direct component. The respective time 
evolution is determined by the fast single ex- 
ponential kinetics (see Figs. [7]and[§]). This 
behavior is related to the small population of 
the charged molecular states. The conclusion is 
that at an off-resonant regime, the appearance 
of a large switch-on current (in comparison to 
its steady state value) could be only related to 
one of the sequential components but not to the 
direct component (compare Figs. [7] and O with 
Fig. [8]). 



B. Resonant regime of charge 
transmission 

The resonant regime of charge transmission 
is formed if the transmission gaps AE + * and 
AE-* become negative. The temporal evolu- 
tion of the respective current components is 
represented in the Figs. [T0]-[T2j The used light 
intensity is identical with the one taken for the 
study of the off-resonant regime. 

1. Single- channel resonant regime 

This regime is achieved if E* > E + + Ep or 
E* + Ep > i.e. if AE+* < or A£L, < 0, 
respectively (cases (c) or (d) of Fig. [5]) . Phys- 
ically, the cases AE + * < 0, AE-* > and 
AE + * > 0, Ai?_* < do not differ from each 
other. Therefore, for the sake of definiteness, 
let us analyze the case (c). As in the previous 
subsection, a situation is considered where the 
fast and the slow kinetic phases are clearly dif- 
fer from each other. Accordingly, the relaxation 
of the current components to their steady state 
values is described by eqs. (l50l) - ([55)l where now 
AE + * < and, thus, 





TIME (ps) 



FIG. 10: Single-channel resonant charge transmis- 
sion with the participation of the charged molecular 
state M+ (the energy gap AE+, is negative). The 
total current is mainly determined by the sequen- 
tial components. 

A comparison of Fig. [10] and Fig. [8] shows 
that the change of AE +te from positive values 
(AE + * =0.1 eV) to negative ones (A£+* = 
—0.1 eV) results in a significant increase of the 
sequential and distant current components (de- 
spite the fact that the width parameters are 
taken even less than those used in the Figs [5] - 
ITD|) . At the same time, the ratio ([52)1 between 
the switch-on and the steady state sequential 
current component is conserved. In the case of 
the resonant regime, one can also introduce the 
ratio defined by the direct current component: 

mr = ijSp/ig ■ (57) 

In line with the expressions (|54|) and ([55]) it 
yields r/dir =^l /^h ■ If one takes the same re- 
lation between the width parameters as it has 
been used in Fig. [8] then rj dlr ~ /T^ ~ 10 
in correspondence with the exact results de- 
picted in Fig. [10] The above given results re- 
fer to a current formation connected with the 
transmission along the channel which is asso- 
ciated with the charged molecular state M + . 
The channel includes two types of transmission 
routs depicted in Fig. [5] the left sequential 
route (M* — > M + — > Mo) and the direct tunnel 
route (M* — > Mo). Analogously, one can con- 
sider the current formation if a charge transmis- 
sion occurs preliminary along the channel re- 
lated to the molecular charged state M_ (right 
sequential route M* — > M_ — > Mo and direct 
tunnel route M* — )• Mo). Recall that the dis- 

(12) (21) 

tant rate constants Ql as well as the Ql 
are defined by both charged molecular states 
M + and M_ . To derive respective analytic ex- 
pressions, one sets K + * = in the eqs. ([40]) . 
([41]) . and P^]) . This results in an analytic 
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TIME (ps) TIME (ps) 



FIG. 11: Single-channel resonant transmission with 
the participation of the charged molecular state 
M-. 



form which follows from eqs. (j43|) and ([55]) if 
one replaces -RT_*( ) an d if*(o)_ by if +H ,( ) and 
-^*(o)+j respectively. Our studies show that in 
this case, the direct current component is com- 
parable with the sequential components (see 
Fig. ITTj) . Since the fast regime of charge trans- 
mission is associated now with the hopping of 
an electron into the HOMO, the maximal value 
of the sequential current component is less than 
that of Fig. GHue to the condition < . 



2. Two-channel resonant regime 

This regime is realized if both transmission 
channels associated with the molecular charged 
states M + and M_ participate in the electron 
transfer process and if the respective transmis- 
sion gaps AEL* and AE + * are negative (cases 
(d) and (e) in Fig. [5]). In this two-channel 
resonant regime, the left and the right electron 
transfer channels represented in Fig. [5] give a 
comparable contribution to the current. The 
time-dependent evolution of the current com- 
ponents is described by the general expressions 
and (HIT) and the set of kinetic equations 
([34]) . Fig. [12] does not show any different be- 
havior among the particular currents which be- 
long to a particular channel. The insert of Fig. 
[T2l demonstrates that the fast part of the time 
evolution completely corresponds to the kinet- 
ics of formation of the excited molecular state 
M* . The slow part reflects the kinetics at which 
a population of the charged molecular states 
M_|_ and M_ varies due to a depopulation of 
the light-induced state M*. Such a depopula- 
tion is negligible during an off-resonant regime 
of charge transmission but becomes pronounced 



FIG. 12: Two-channel resonant charge transmis- 
sion. There is no basic difference between the 
behavior of the current components formed by a 
single-channel of charge transmission (except a cer- 
tain increase of the direct component, cf. also Fig. 
1 10p . The evolution of the current components fol- 
low the time-dependent behavior of the molecular 
probabilities (see insert). 



in the resonance regime. 



VI. CONCLUSIONS 

We put forward a detailed study on the time- 
dependent behavior of light-induced transient 
currents in molecular junctions, like in a molec- 
ular diode l-M-2. A nonequilibrium set of ki- 
netic equations has been derived for the molec- 
ular states which participate in the current for- 
mation. Those are the neutral molecular states 
Mo and M* and the two charged molecular 
states M + and M_. It could be shown that 
an interelectrode electron transfer 12 <=± l + 2 
takes place along the channels associated with 
the charged molecular states M + and M_ , see 
Fig. [5] These states participate in a light- 
induced interelectrode electron transfer in the 
absence of an applied voltage either as real in- 
termediate states (forming the sequential trans- 
mission route) or as virtual intermediate states 
(forming the direct transmission route). The 
sequential route includes the hopping of an elec- 
tron between the molecule and the adjacent 
electrodes, being thus responsible for a molec- 
ular charging. The formation of the respec- 
tive direct current component is accompanied 
by the transition of the molecule from its pho- 
toexcited state M* to its ground-state M . The 
time-dependent behavior of the total photocur- 
rent is governed by the molecular populations 
P(a;t), cf. eqs. (0DJ and (gTJ). The latter 
evolve in line with the kinetic rate equations 



16 



(EH)- 

The relative efficiency of each route is de- 
termined by the actual value of the transmis- 
sion gaps AE a *, (a = +,—): If AE a * is pos- 
itive, then an electron transmission along the 
M Q -channel proceeds in an off-resonant regime 
whereas for a negative AE a * the current is 
formed in the resonant regime. In this regime 
the sequential and the distant current compo- 
nents significantly exceed the same components 
formed at the off-resonant charge transmission 
(see in this context the Figs. [TU] and [11] and 
compare them with the Figs. [7] and [8]). In 
the framework of a HOMO-LUMO model, the 
direction of the light-induced electron current 
is determined by the factor T^r^ - T^r^ 
which reflects the difference between interelec- 
trode electron flow 12 — > l + 2~ and 1 2 + «— 
12. 

If the difference between the width param- 

(r) 

eters Ij- becomes large, a characteristic ki- 
netic effect appears for the transient photo cur- 
rents. In this case, those currents do exceed 
their steady state value to a large amount. The 
physical origin of this effect is related to the 
fact that the photoexcitation of the molecule 
(Mo — > M* process) initiates a molecular charg- 
ing (see the schemes in Fig. @] and Fig. [5]). 
Charging is caused by the transition of an 
electron from the molecule to each electrode 
(Af» — > M + process) or from the electrodes to 
the molecule (M* — > M_ process). Such a light- 
induced electron motion forms the fast initial 
part of the electron transmission which can be 



seen in the time-dependent behavior of the se- 
quential current components Ise\(t) and lieq(t). 
Molecular recharging is characterized by con- 
tact (hopping) rate constants. If the character- 
istic recharging time is much less than the char- 
acteristic time r st of the steady state formation, 
then the maximal value of the transient pho- 
tocurrent may become rather large compared 
to its steady state value (see, for instance, Fig. 
[inland Fig. HU). 

Experimental studies of transient photocur- 
rents allow one to clarify the details of contact 
(molecule-electrode) and distant (electrode- 
electrode) electron transfer processes in molec- 
ular junctions. We have shown that the ef- 
fective formation of the photocurrent becomes 
possible if the photon energy ftw — E* — Eq ex- 
ceeds the energy gaps AE_ or/and AE +0 , cf. 
Fig. [2]). This corresponds to an energy level 
arrangement as shown in the schemes (c)-(f) of 
Fig. El 
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